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^ ' Abstract 

o ■ 

The six-quark instanton induced 't Hooft interaction is considered in combination with the Nambu- 
Jona-Lasinio (NJL) type Uii'i) x t/ij(3) chiral symmetric Lagrangian. We discuss the bosonization 
^ , of this multi-quark interaction, taking the Ua{1) breaking as a perturbation. We discuss its relation 
^ I with the usual approach. 

> , 

Q\ . Introduction 
0\ • 
(N ■ 

■ In low-energy QCD phenomena one considers basic powerful classification schemes, as symmetries 
O ■ and symmetry breaking patterns, large Nc counting rules PP- |1] and scales at which some type 
Q . of interactions set in j3] . A well known example is rj' meson physics, where all these concepts are 
"Th '. needed to explain its large mass m^/ ~ IGeV. If the singlet SU{3) axial current of QCD, 
Ph! were conserved as m^urr ~^ 0, the rj' would occur as the ninth pseudoscalar Goldstone. As the 
o il global f/L(3) X Ur{3) chiral symmetry of QCD is broken by the Ua{^) Adler-Bell-Jackiw anomaly 
^ ■ jH] , it was realized through the study of instantons [7j-[n] that effective 2Nf quark vertices arise, 
^ . the 't Hooft interactions , responsible for the non-conservation of the singlet axial current and 
• ^ ! absence of the related Goldstone boson. Typically, the mass squared of an approximate Goldstone 
^ I boson is linear in the symmetry breaking parameter. The anomaly is of order 1/Nc, one expects 
c3 ] m^/ to be of order (l/Nc). Although the Ul{S) x Ur{3) symmetry is recovered in the limit Nc — > oo, 
the empirical mass of the rj' is abnormaly large to conform alone with the large Nc counting rules, 
which is Witten's puzzle ^U]. The solution was given in j3], where with the help of QCD sum 
rules it was shown that a large scale emerges in the 0~,0"'" gluonic channels, compared to which 
the r]' mass is small, and characterizing the breaking of asymptotic feedom in these channels. 
Our plan consists in addressing the rj' physics and lowest lying (0~, 0"*^) meson nonets in a many- 
fermion vertices model without explicit gluon degrees of freedom. In this way one takes instead 
the quark structure of the mesons explicitly into account, and may obtain relevant information 
about the vacuum structure of QCD, albeit in an effective field theoretical approach. We use the 
successful Nambu - Jona-Lasinio model ^I] extended to the U{3)l x U{3)r chiral symmetry of 
massless QCD and the six quark 't Hooft interactions for the symmetry breaking. We review 

some known bosonization schemes with functional integral methods and show their restrictions. 
Within these schemes the leading order t]' mass is obtained close to its empirical value and a good 
description of the remaining members of the lowest pseudoscalar nonet and related weak decay 
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constants and quark condensates. The i]' puzzle is also present in this case. Finally, using the 
spectral decomposition approach, we indicate how to obtain systematically all f/A(l) corrections 
in a perturbative way and its relation to the previous approach. 



The model Lagrangian 

The starting Lagrangian is 

C = q{i-i^d^ - rh)q + Cnjl + Cq, 
with NJL four-fermion vertices of the scalar and pseudoscalar types 

G r 



C-NJL — — 



{(iKq) + {qiibKqy (2) 

and the six-quark 't Hooft interaction 

Cg = K (det qPLQ + det qPRq), (3) 

Here the positive coupling G, [G] = GeV~^ has order G ~ l/N^, the negative coupling k, [k] = 
GeV~^ with large Nc asymptotics k ~ . The Pl^r = (1 =f 75)/2 are projectors and the 

determinant is over flavor indices. At large Nc-, Cnjl dominates over Cq. 

This Lagrangian has been studied at mean field level in terms of quark degrees of freedom in ^21 
and in a bosonized form in ^2] and has been widely used since then ^3]. . 

Bosonization 

We use functional integral methods to bosonize the multi-quark Lagrangian, for which the vacuum 
persistance amplitude is 

Z = j VqVq exp (^i J d^x . (4) 
The fermion vertices of the original Lagrangian are linearized with help of the functional identity 
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with 



/ YiVSaVpaSiSa - q\aq)KPa " qilbKq) 
a 

X exp |z j d^x [aaSa - q\aq) + (t>a{Pa - qil^Kq)] 

Z = j '[[VaaV(f)aVqVqexp (^i J d'^xCg{q,q,a,(p) 

^ J Yl'I^Sa'Dpa exp (^i J d^xCr{(y,(t>,s,p)^ (5) 

Cq{q, q, a, (f)) = q^i-f^'d^ - a - i-i^(t))q, (6) 
G 

Cr{a,(j),S,p) = — [sI+pI^ + SaicTa- rha) +Pa(f)a 

+ — Kbc Sa {SbSc - ^PbPc) ■ (7) 



The totally symmetric constants Aabc are defined through the flavor determinant 

det U = AabcUaUbU, , (8) 

where U = UaXa- The scalar fields cTq, a = 0, 3, 8 have non-zero vacuum expectation values related 
to spontaneous breaking of global chiral symmetry, which give rise to the constituent masses 
of the quarks. New scalar fields with zero vacuum expectation values (0|(Ta|0) = must be defined 
through the shift cr^ — s> + rria. We seek the final bosonized Lagrangian expressed in terms of 
tree level mesonic fields cxa and 0a and must integrate out all remaining fields. 



Methods used for the functional integration 



The Gaussian integration over the quarks is done using the generalized heat kernel technique 
developed in ^21- It accounts for a chiral covariant treatment of the one- loop determinant of the 
Dirac operator with an arbitrary non-degenerate quark mass matrix at each order of an expansion 
with generalized Seeley-DeWitt coefficients. 
The integral over the auxiliary bosonic fields Sa,Pa 

I[aa,(j)a] = j YlVSaVpa CXp (^i J rf^X £,.((7, 0, S, p)^ . (9) 

is done perturbatively. For that let us devide the lagrangian Cr in two parts. The free part, Cq, 
given by 

G 

Coi(^A,S,p) = — [sI+pI^ + Sa(Ta+Pa^a- (10) 

The 't Hooft interaction is considered as perturbation £/ 

jCl{Sa,Pa) = :^AabcSa{SbSc-3PbPc) ■ (H) 

For simplicity, a one-dimensional field theory approximation for the integral is considered. 
The extension to the real case is straightforward. We put our system in the interval of size L, i.e. 
—L/2 < X < L/2, assuming the limit L oo in the end of our calculations. 
The Fourier decomposition of the fields fa{x) = {sayPayO'ay'Pa} inside the interval is 

fa{x) = E fn exp (.2vr^) . (12) 

71= — OO ^ ^ 

This corresponds to the periodic boundary conditions /„(— L/2) = fa{L/2). Then we have 

f fL/2 

J dxCr{x) ^ j dx Cr{x) = L Cr , (13) 

where Cr = '^o + '^i, and 

+ 00 r 

^ ' ^ (14) 

1^1 = ^^<^bc E - ^pLpI) Sn+m+k,0- (15) 



+00 r 

^0= E ^ (<^-n + P>-„) + + CP-, 

n=— 00 '- 
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The functional integral can be understood as the product of integrals over the Fourier coefficients 

+00 +00 

/ Vsa{x)Vp,{x)^ll f dsldpl. (16) 

—00 —00 

Thus, in the perturbative approximation we have 

+00 

/[(T,0]~exp{zL£,(5:,P,")}n / dsldple^v{iLCo) (17) 

n •' 
—00 

where 



It is convinient to use the normalized functional /7v[cr, 0] defined as 

/.K*l^||. (18) 

The Gaussian functional integrals can be simply evaluated, for instance, one has 

[ '^'^ ( G \1 

-^M = n \ C^SnCXpjzL ^ — + Cr„,S_„j I 



iL 



-00 



/[0]exp -— ^n0^n\- (20) 



Thus the functional integration yields 



/;v[a,0] ~ exp{zL£,(^:,P:)} 



X ^^P|~^ 5Z (o-nO'-n + </>n</'-n)| • (21) 

One should calculate partial derivatives to obtain the effective action Fe// 

J^[a,0] ~A[a,0]e*^^//. (22) 

We have 

Fe// = i In 0] + Fq - z In (l + e"^^" (e^^^ - l) e^^") , (23) 

where F/ = Lti{S^, P^), and A[cr, 0] is fixed by the requirement that the effective action Fg^j is 
real. Here Fq represents the leading order action. 

T OO 

rO = -^ E (<^-n + Cr-n)- (24) 
^ n=— OO 

The logarithm in eg. ()23j) is a source of f/(l) breaking corrections which arise as a series in powers 
of the partial derivatives. Expanding 

OO 

5 = e-'^» (e'^^ - l) = ^ (25) 



n=l 
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one can determine 



5i 



-iL 
32G3 



(26) 



n,m,k 



r _ . 9L ^ 

^2 — 2 * 64(2(7)^ ^abc^abc ^m+k+fh+kfi 



m,k,m,k 



X 



(8G2)2 ^ 



+ 



3z 



El 



(27) 



32X^3 

and so forth. The factor A[a, 0] must be also expanded 

A[a, (f)] = l + K(3i + K^f32 + 0(ft:^), 
and, up to the considered order in n, we have 

A = 0, P2 = E (<^-n + K<pQ E 1 • 

Therefore the perturbative action systematically obtains U{1) breaking corrections in k", which 
we collect in the corresponding part of the action, r„. 



where Tq is given by eq. (j24j) . and 



32^3 



+ 00 

^eff = E , 
n=0 



^abc E (^ni^^m'^k - ^<PL<Pk)^n+m+k,0 



(28) 



n,m,k 



64(2G)! 



A-abcA-abc E ^m+k+ni+k,{ 



m,k,m,k 



[^rh^k 



k - ^'k<Pl 



(29) 



Here an unessential constant has been omitted in Taking the infinite - volume limit L oo, 
one finally obtains 



+ 00 



K 



326-3 



64(2G) 



2G I i^aix) + (l)lix)) , 

— oo 

■oo 

— OO 

+ 00 

AabcAabc dx { 4ab{x)ai{x)(j)c{x)(j)c{x) 
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+ [(Tb{x)acix) - (f)bix)(t>c{x)] [ai{x)ac{x) - (f)i{x)(f)c{x)] } 
and so forth. The infinite volume hmit for the coefficient (32 yields 

rL/2 1 +00 

L/2 



1 fL/2 I +°° 

^2 = 7^7^ / dx {aaix)aa{x) + (f)aix)(f)a{x)) y I] 1 • 
8G^ V= J~L/2 L 



The integral has a smooth L oo limit, but contains a local factor, which can be understood as 
5-function singularity, 5(0), see also ^H]; To show this one needs the following result 

^ exp U-^n \ =L ^(3^ - = L^i^)- (30) 

?i=— 00 \ / 

On the last step we took into account that in the considered problem all x-dependent functions 
are integrated only in the interval — L/2 < x < L/2 and, therefore, only the term with n = 
contributes. The infinite value 

1 +00 

7 E 1 = ^(0) (31) 

represents the density of Fourier harmonics in the interval. It must be regularized by cutting an 
upper part of the spectrum, e.g. 

1 ^ 27V + 1 
^(O)reg = 7 E 1 = (32) 

^ n=-N ^ 

where is large enough. One can fix ^ 1 without any relation to the size of the box. In this 
case the limit L ^ 00 leads to the vanishing of the (5-function and to = 0. Alternatively, one 
can relate A^ with L by introducing the momentum space cutoff A^;: N{L) = LA£;/(47r) ^ 1. 
Unlike L the cutoff A^; has an obvious physical meaning giving the scale of momenta relevant for 
the problem. Indeed, the nth harmonic has a momentum p„ = 2im/ L. The size of the considered 
box is the difference — P-n = 4:nN/L = Ae- This scale cannot be eliminated by taking the 
limit L — s> 00. One has instead 

1 

27r L 

where only the second term does not contribute in the limit L — > cxd. By introducing the cutoff 
A^;, we suppose that the density of Fourier harmonics has a finite value which can be fixed 
phenomenologically. In this case ^ 0. It is a pure physical matter of argument to choose 
among these two alternatives. 

The Reinhardt and Alkofer estimate ^3] takes the integral (jOJ at one stationary phase point 

without any integration over the auxiliary variables Sa,Pa- 

The contributions at the stationary phase condition are obtained as 

^r(rst) = § tr {U,tUl) + ^ tr {WUl + W^U,t) 
1 1 

= haCa + -hab(^a(rb + l^h'^ab <P a(pb + • • • • (34) 



'^(O)rcg = 7^ + T ' (33) 



fulfiUing 



GUa + Wa + ^AabcUlUl = (35) 
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with covariant combinations 

Wa = aa + Aa- i(t)ay 

and definitions 

Aa = - rha, Ua = Sa - ipa- 

The field Ust represents the exact solution of condition (jH^ . obtained by expanding Sa,Pa in parity 
even and odd combinations respectively, of increasing powers of bosonic fields (pa, c"a 

s^st = ha + h'-^J ab + h'-^labac + h'^abcMc 

+ h^abcd^b(^c(Td + h^abcd^h^c(pd + ■■■ 
pit = h'ab (Pb + h''al(l)b(Tc + h^abcd(^b(Tc(pd 

+ h'^ldMccPd + .... (36) 

The coefficients /i^^'' depend explicitly on the quark masses and coupling constants G, k and are 
fixed by a series of coupled equations following from (jH^j) and obtained by equating to zero the 
factors before independent combinations of mesonic fields. Due to recurrency of the considered 
equations all coefficients are determined once the first one, ha, has been obtained. Note that the 
stationary phase condition has several solutions [T7]. If one wants to apply the stationary 
phase (SP) method consistently, one must take into account the contributions of all critical points 
(even at leading order). The 1/Nc counting indicates that the last term in the stationary phase 
eq. (j35| is l/Nc suppressed. This argument can be used for the systematic l/Nc expansion of the 
effective action, what actually has been done in our perturbative approach. 

We have found strong indications that the estimate PSI corresponds to a resummation of the per- 
turbative series in k, collecting all tree (non-singular) contributions related with the auxiliary fields 
Sa,Pa of the perturbative approach. Further terms of the ressummed series yield systematically 
"loop" -corrections of these fields pOj . 

Results and Discussion 

We refer to [TH] for a detailed derivation and calculation of the leading order low lying pseudoscalar 
and scalar meson spectra and related observables. Further sets are also presented in (TH] . Already 
in its simplest form, the Lagrangian yields overall quite good fits, even the scalars are reasonable. 
They have been obtained with the following approximations. 

a) Related with the evaluation of the Gausiian functional integral over quarks: knowing that the 
lack of confinement in the NJL model introduces serious difficulties with the crossing of non- 
physical thresholds associated with the production of free quark - antiquark pairs, which one 
may encounter by formally continuing the full Euclidean action to Minkowski space, the heat 
kernel series has been truncated at second order. We admittedly deviate from the original NJL 
Lagrangian, however in a way which relies heavily on its symmetries and asymptotic dynamics, 
which are fully taken into account. 

b) Related with the functional integral calculation over the auxiliary variables Sa,Pa, of cubic 
order: we have used the SPA with one critical point. We know however that there exist more 
critical points which should be taken into account even in leading order in SPA. 

Finally we have suggested a different way of bosonization, considering the t Hooft determinant 
as a perturbation to the four-quark intercation. One can relate the perturbative expansion to 
the SPA with one critical point through a resummation. This is discussed in full detail in |2Uj . 
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Corrections to the leading order can be obtained in a systematic way. Due to the local character 
of the interactions these corrections involve a divergent density of harmonics. We discuss how a 
cutoff can be introduced, and together with the ultraviolet cutoff related with the quark loops, we 
obtain an effective bosonized Lagrangian. Calculation of the related mass spectra is in progress 

m 
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